F=rand(3,3);
[U,D,V] = svd(F,0);
FB = U*diag([D(1,1) D(1,1) 0])*V';

FBE=FB+rand(3,3);




  
    [V,D] = eig(FB);
    B{j}=V(:,1:2);






for i=1:numTries
    focal=fstart+(randn()*fstd);
    ar=(1+rand()*arstd);
    skew=skewstd*randn();
    xguess=xstart+(randn()*xstd);
    yguess=ystart+(randn()*ystd);
    Kguess = [focal   skew        xguess;   0      focal*ar   yguess;   0            0             1  ];
    
    
    for j=1:numFs
        G=(Kguess')*F{1,j}*Kguess;
        [V,~] = eig(G'*G);
        vv= subspace(V(:,1:2),B{j});
        Ed(i,j)=vv(1,1);
        
    end
    
end